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Abstract. The production of tj-mesons in the pp ppui reaction has been investigated with the COSY- 
ANKE spectrometer for excess energies of 60 and 92 MeV by detecting the two final protons and recon- 
structing their missing mass. The large multipion background was subtracted using an event-by-event 
transformation of the proton momenta between the two energies. Differential distributions and total cross 
sections were obtained after careful studies of possible systematic uncertainties in the overall ANKE ac- 
ceptance. The results are compared with the predictions of theoretical models. Combined with data on the 
meson, a more refined estimate is made of the Okubo-Zweig-Iizuka rule violation in the 4>/^ production 
ratio. 

PACS. 13.75.-n Hadron-induced low- and intermediate-energy reactions and scattering (energy < 1 GeV) 
- 13.60.Le Meson production - 14.40.Cs Other mesons with S = C = 0, mass < 2.5 GeV 



1 Introduction 

The nucleon-nucleon system is strongly coupled to chan- 
nels that contain one or multiple mesons. The production 
of these in iVA^-collisions, preferably near their respec- 
tive thresholds, will test and constrain theoretical models 
since, in these cases, only few partial waves contribute. 
The production of pions and heavier pseudoscalar mesons 
including kaons has been systematically studied in proton- 
proton interactions in the near-threshold region [1]. Now 
vector mesons, such as the p and the w, are significant con- 
tributors to the NN force at short distances but for these 
particles much less information is available on their pro- 
duction mechanism due to the lack of experimental data. 
The relative production of the isoscalar uj and mesons 
has recently received renewed interest in connection with 
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the so-called Okubo-Zweig-Iizuka (OZI) rule [2]. This 
ratio may give information on the admixture of strange 
quarks in the nucleon only if the production mechanism 
is known to be the same for both mesons [3-5]. This can 
be best controlled through near-threshold measurements 
but, as yet, the experimental data set is rather limited. 

In the m case, the pp ppuo total cross section was 
measured for excess energies Q < 30 MeV at SATURNE 
by detecting the two protons in the SPESIII magnetic 
spectrometer and identifying the meson by the missing- 
mass method [6]. In a more exclusive experiment, the 
DISTO collaboration also measured oj production but at 
much higher energy {Q « 320 MeV) [7] . The only data be- 
tween were taken at COSY-TOF, a non-magnetic time- 
of-flight spectrometer, where the reaction was studied at 
Q = 92 and 173 MeV [8]. This latter energy is in fact 
the lowest for which differential distributions are avail- 
able. Extra data are therefore required to fill the gap and 



2 



S. Barsov et.al: Study of w-meson production in pp collisions at ANKE 



confirm the energy dependence. For this purpose we here 
present co results obtained using the COSY- ANKE facility 
at 60 and 92MeV. It is relevant within the OZI context 
to note that the DISTO group measured d> production at 
82MeV, and we have recently published similar data at 
18.5, 34.5, and 75.9 MeV [9]. 

The experimental facility at our disposal is described 
in section 2, where particle identification and efficiency 
determination, as well as the extraction of the luminos- 
ity, are also discussed. Since the w-meson production is 
identified from the missing mass of two final protons, one 
of the major problems is a large background under the lu 
peak arising from the production of two or more pions. 
We use the same kinematic transformation as proposed 
in the SPESIII analysis to estimate this multipion back- 
ground by employing data obtained at the other energy. 
This is the crucial element of the pp ppuj event selec- 
tion, which is treated in detail in section 3. The acceptance 
of the two protons from the pp — > jypuj reaction in ANKE 
was far from complete and so a thorough investigation 
of various one dimensional differential distributions was 
undertaken and compared with the results of simulations 
obtained on the basis of various model assumptions to be 
found in section 4. The values of the total cross sections 
given in section 5 are compared with the results of several 
theoretical approaches. Also to be found there is a discus- 
sion of how our data influence the extraction of the OZI 
ratio from near-threshold oj / (j) production. Our conclu- 
sions and thoughts for future investigations are presented 
in section 6. 



2 Experimental set-up and raw data analysis 

The measurement has been carried out at the ANKE in- 
stallation [10] , using the internal hydrogen cluster-jet tar- 
get [11] and COSY circulating proton beam. To reduce 
systematic effects, the beam momentum was changed ev- 
ery 10 minutes between 2.85GeV/c and 2.95GeV/c, cor- 
responding to excess energies of 60 MeV and 92 MeV, re- 
spectively. Charged particles resulting from beam -target 
interactions, and going forward with laboratory polar an- 
gles up to 20°, were deflected by the central dipole D2 
onto the ANKE detector systems. In the present investiga- 
tion, one final proton was detected in the ANKE forward 
dc;te(;tor (FD), passing through a set of three multiwire 
proportional chambers (MWPC) and two layers of scin- 
tillators placed downstream of D2 and close to the beam 
pipe. The other proton was detected in coincidence in the 
positive charge detector system (PD), positioned to the 
right of D2. This system includes one layer of 23 thin start 
scintillators (SA) placed close to the exit window of the 
D2 vacuum chamber, two MWPCs, and one layer of six 
large stop scintillators (SW) within the ANKE Side Wall. 
Only that part of the PD which has acceptance for the 
pp — > ppu) reaction was used in this experiment. 

In order to suppress the very large count rate from 
proton-proton elastic scattering, the coincidence of sig- 
nals from the FD and SW scintillators was chosen as the 
main trigger. However, 1% of the count rate from the FD 



scintillators alone was added to the trigger stream so that 
the data could be normalised. The momenta of these par- 
ticles were reconstructed and pp elastic scattering identi- 
fied from the clear peak in the missing mass distribution 
at the nominal proton mass. The data were corrected for 
the efficiency of the FD MWPCs on an event-by-event 
DSrSlSj clS described below. The number of protons elas- 
tically scattered between laboratory angles of 6.6° and 
8.8° was determined from a fit of the peak by a Gaussian, 
with a polynomial background. The contribution of the 
background under the peak was found to be (1.8 ± 0.8)%. 
Using predictions for the differential cross section from 
the SAID database [12], the same average luminosity of 
7.5 X 10^°cm~^s~^ was found for the two incident mo- 
menta, the values being constant over the angular range 
to better than 2%. The corresponding integral luminosi- 
ties were 0.280 pb"^ and 0.273 pb'^ at 2.85GeV/c and 
2.95GeV/c, respectively. The overall systematic error of 
about 5% was estimated taking into account the uncer- 
tainty in the SAID calibration standard (~ 3%), the ef- 
ficiency determination (~ 1%), and the 0.1° precision in 
the determination of the scattering angle 3%). 
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Fig. 1. The time-of-flight distributions of particles detected 
in one of the SA and one of the SW scintillators, obtained 
before (dashed) and after (solid line) the determination of their 
momenta. The reduced number of events in the latter case is 
the result of the cfRcioncics of the MWPCs and the track- 
finding algorithm. The left and right peaks arise, respectively, 
from pions and protons emitted from the target. 



As the pp — > ppuj reaction has to be identified from 
the missing mass with respect to the two protons, these 

were distinguished from other particles, mostly pions, us- 
ing the time-of-flight (TOF) technique and momentum 
determination. 

Although the SW scintillators arc capable of deliver- 
ing a momentum acceptance of the PD up to 1.6GeV/c, 
in the actual configuration of the equipment the momen- 
tum range had to be limited to be below l.lGeV/c. Un- 
der these conditions, the time of flight between SA and 
SW scintillators alone provided sufficient particle identifi- 
cation. As an example of this, Fig. 1 represents the TOF 
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spectrum corresponding to the worst case of pion-proton 
separation. 

The total efficiency of the momentum reconstruction 

in the PD could be determined for each of the SA-SW 
combinations because both of the MWPCs are placed be- 
tween these scintillators. The numbers of protons before 
and after the reconstruction of their momenta were ex- 
tracted from the corresponding TOF spectrum, which is 
similar to the one presented in Fig. 1. The accidental back- 
ground was significantly suppressed by the requirement to 
have some particle track reconstructed in the FD. The ef- 
ficiency was found to decrease smoothly from 92% to 86% 
with increasing particle momentum and track inclination 
but to remain constant in the vertical direction for all 
SA-SW combinations involved in the later analysis. 

' J 




FD momentum [GeV/c] 



Fig. 2. The TOF difference between particles detected in the 
FD and protons in SW scintillators versus the momentum of 
particles in the FD. Protons detected in the FD produce the 
upper curved band whereas the TOF of pions is almost inde- 
pendent of the momentum. 



In contrast to the PD, the two layers of scintillators 
used in the FD are placed very close to each other and, 
in order to identify protons, the TOF difference between 
the SW and FD scintillators had to be combined with the 
measurement of the momentum of particles detected in 
the FD, as shown in Fig. 2. Coincident protons were se- 
lected in the PD as described above and their momenta 
reconstructed to make a correction for the different trajec- 
tory lengths. In the momentum range of interest, which is 
below 1.8GeV/c, the proton and pion bands are well sep- 
arated except in the vicinity of 1.75 GeV/c. Here the num- 
bers of pions are strongly enhanced by the proton neutron 
final state interaction in the pp — > pmr~^ reaction. How- 
ever, if misidentified as protons, such pions will produce a 
peak in the pp missing-mass distribution at 0.65GeV/c^ 
and 0.68GeV/c2 for 2.85GeV/c and 2.95GeV/c, respec- 
tively. This allows the possibility to control the proton 
selection also in the overlap region of Fig. 2. The total 



loss of protons as a result of the selection was estimated 
to be (7 ± 2)%. 

To determine the efficiency of a single FD MWPC, the 
trajectory of a particle has been found using only infor- 
mation from the other two MWPCs. Having calculated 
the point at which the trajectory intersects the sensitive 
plane of the MWPC under investigation, the presence of a 
valid neighbouring cluster has been checked to satisfy all 
requirements applied in the track-finding algorithm [13]. 
This procedure resulted in an efficiency map for each of 
the sensitive planes, which was subsequently used to intro- 
duce the efficiency correction on an event by event basis. 
Only one sensitive plane among the seven involved in the 
FD track reconstruction was found to have an inhomoge- 
neous distribution of efficiency in approximately 30% of 
its sensitive area, with an average value of 88%. All the 
others had homogeneous maps with average efficiencies 
above 96%. 

The influence of the total efficiency corrections on the 

missing-mass distributions has been checked using the ra- 
tio of corrected to uncorrected spectra, which was found 
to be constant to within experimental errors. The average 
efficiency was about 75%. 

3 Selection of pp ppw events 

Histograms representing the overall missing-mass (mx) 
distributions from the pp ppX reaction at the two beam 
momenta are shown in Fig. 3. 




mx(p,p) [GeV/c^] 

Fig. 3. The missing mass distributions of the pp ppX reac- 
tion measured (a) at 2.85 GcV/c and (b) at 2.95 GcV/c. There 
is clear evidence for peaks corresponding to the production of 
both the ri and lj (indicated by arrows) but sitting on a large 
multipion background. 



Both spectra show similar features. There is a sharp 
rise from the kinematic limit on the right due to the mul- 
tiparticle phase space, whereas the more gentle fall to the 
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left is mainly a reflection of the ANKE acceptance. Sit- 
ting on top of this multipion background there are small 
peaks due to the production of the r]- and w-mesons. Ev- 
idence for the pp ppn'^ reaction was also found below 
0.4GeV/c^, but the pion peak has a large width because 
the effects of momentum uncertainties are significantly 
amplified in this region of large excess energies. 

The accuracy of the missing-mass reconstruction was 
verified using the pp pp-q and the pp pmr^ reactions, 
with the results being presented in Table 1. In the pion 

Table 1. Missing masses of various reactions in MeV/c^ for 
the two beam momenta. 



case, the neutron mass was reconstructed in both possible 
variants, cither when the proton was found in the FD and 
the pion in the PD or vice versa. Given the good agree- 
ment achieved in these cases, one should expect to find the 
uj peak at its nominal position. Howcn^cr, in the vicinity 
of the LU mass, the large background from multipion pro- 
duction is strongly peaked, due to the specific acceptance 
of the ANKE magnetic spectrometer, and the challenge 
is to extract the number of pp ppui events. To select 
the LU signal from Fig. 3, a kincmatical transformation 
of the experimental missing-mass distribution from one 
beam momentum to the other was applied, as illustrated 
in Fig. 4. This approach was proposed and successfully 
used in the study of w production at lower excess energies 
with the SPESIII spectrometer [6]. It avoids the need to 
construct specific models for a background that is unlikely 
to follow phase space. The scheme involves an event-by- 
cvent Lorentz transformation of the momenta of both final 
protons from the laboratory system where they were actu- 
ally measured to the laboratory system at another beam 
momentum. If the phase space of the multipion produc- 
tion processes is much larger than the experimental accep- 
tance, the shape of the missing mass distribution related 
to this background is expected to be fixed by the accep- 
tance and to remain the same in both systems. On the 
other hand, a peak corresponding to the production of a 
single meson is shifted by approximately the difference in 
centre-of-mass energies, as seen clearly for r] production 
in Fig. 4. The solid histogram in the top panel of this 
figure represents the missing-mass distribution actually 
measured at 2.85GeV/c (see Fig. 3a). The dashed his- 
togram represents the result of such a transformation of 
the missing-mass distribution obtained at the 2.95 GeV/c 
beam momentum (see Fig. 3b) to the 2.85GeV/c. The 
transformed distribution is normalised by the ratio of the 
sums of events in the 0.6-0.7 GeV/c^ range, which gives a 
factor of 1.054 ±0.006. 
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Fig. 4. Upper frame: The missing mass distribution of the 
PP ~* PP^ reaction measured at 2.85 GeV/c (solid line, Fig. 3a) 
is shown together with the analogous distribution obtained at 
2.95GeV/c (Fig. 3b) and kincmatically transformed to the 
2.85GeV/c (dashed line), as described in the text. Lower 
frame: The bin-by-bin difference of these two distributions is 
represented by points. The solid lines show simulations of the 
pp — » ppr] and pp ppuj reactions convoluted with acceptance. 
The widths of the peaks, which are strongly influenced by ef- 
fects of momentum reconstruction, are very well reproduced. 



The difference between the actual and transformed dis- 
tributions, shown by points in the lower frame of Fig. 4, 
demonstrates clear positive (negative) peaks correspond- 
ing to r] and ui production in the actual (transformed) 
cases. These peaks are in a good agreement with the sim- 
ulations of these reactions (solid line). The missing-mass 
distributions simulated at 2.95 GeV/c were transformed to 
2.85 GeV/c and subtracted in the same way as the experi- 
mental one. The r] and w experimental widths arise mainly 
from the uncertainties in the momentum reconstruction. 
Finally, the difference spectrum for mx > 0.65MeV/c^ 
was fitted by two Gaussian functions in order to extract 
the number of events in the lu peak, its statistical uncer- 
tainty, and the position of the peak at Q = 60 MeV (see 
Table 1). 

The same procedure was applied to extract the total 
number of w-mesons detected at Q = 92 MeV. In this 
case the missing-mass distribution measured at a beam 
momentum of 2.85 GeV/c was transformed to 2.95 GeV/c. 



mx 2.85 GeV/c 2.95 GeV /c 

mn{n'^U'PPo) 939.3 ± 0.1 939.6 ±0.1 

m„(ppD, TT+o) 939.2 ± 0.2 939.8 ± 0.2 

mr,{ppD,pFD) 547.8 ±0.1 547.4 ±0.1 

mu^(pPD,PFD) 782.3 ±0.4 782.7 ± 0.4 
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However, the good agreement found for the shapes of 
the background spectra in the region between the r] and 
LO peaks does not necessarily prove that the same is true 
also underneath the u peaks. In principle, up to six pions 
might be produced at our energies, but only those with 
two and three pions in the final state have phase spaces 
much larger than the ANKE acceptance. Although the 
total cross sections for the production of four or more 
pions are at least an order of magnitude smaller, they 
have larger acceptances in ANKE and populate missing 
masses mostly near the kinematical limit. 

The direct way of verifying the shape of the back- 
ground under the oj peak, by interpolating between two 
transformed distributions measured above and below the 
given excess energy, is not possible here since data were 
taken at only two energies. An attempt to describe the 
multipion background with the help of phase space sim- 
ulations, adjusting the relative contribution of reactions 
with different number of pions, was also undertaken. How- 
ever, no set of parameters could simiiltancously provide 
an appropriate reproduction of the background shape in 
different parts of the acceptance, probably due to the ne- 
glect of all resonance effects. We have therefore used the 
following procedure to evaluate the possible variation in 
the shape near the kinematical limit in the missing mass. 




mx(p,p) [GeV/c^] 

Fig. 5. The missing-mass distribution of the pp — * ppX re- 
action obtained at 2.85GeV/c within a restricted part of the 
total acceptance (closed circles) is shown together with the 
corresponding distribution kinematically transformed from the 
2.95 GeV/c to 2.85 GeV/c (crosses). The thick solid line repre- 
sents a fit to the distribution at 2.85 GeV/c while the thin one 
shows the background function under the ui peak found in this 
case. The same fit procedure was applied to the transformed 
distribution and the resulting background function shown by 
the dashed line. 



Due to the characteristics of the ANKE installation, 
the shape of background varies across the total spectrom- 
eter acceptance. Using this feature, an acceptance range 
was found where the shape of background could be di- 



rectly established from the fit procedure. The results from 
this are shown in Fig. 5. Both the measured and the trans- 
formed distributions were then separately fitted, using the 
same function to describe the background, together with 
a Gaussian function for the oj peak. The transformed dis- 
tribution in Fig. 5 was scaled by 1.086 ± 0.011, i.e. the 
ratio of the integrals under the background functions ob- 
tained in this way. The ratio of the background functions is 
not quite constant over the range but, apart from the 5—7 
MeV/c^ interval next to the kinematical limit, the devia- 
tion from unity is rather small. The difference between the 
background functions observed near the maximum of the 
distribution does not exceed 5%, changing sign in vicin- 
ity of 0.75GeV/c^, and staying at the 2% level at smaller 
masses. Since the accuracy of the fit is similar to these de- 
viations, a variation of 3% has been considered as a mea- 
sure of the systcmatics when the number of events under 
the LU peak could be determined only from the difference 
of experimental distributions at different energies. 

The total numbers of tj-mesons deduced from fits to 
the distribution at 2.85 GeV/c shown in the lower frame of 
Fig. 4, and the analogous one at 2.95 GeV/c, were found 
to be: 

5560 ± 540(stat.) ± 660(syst.) at Q = 60MeV, 
5360 ± 480(stat.) ± 720(syst.) at Q = 92MeV. 
The large statistical errors are the result of the subtraction 
procedure as well as the need to describe two peaks of 
opposite sign situated close to each other. 

4 Analysis of differential distributions 

The restriction on the particle momenta registered in the 
PD leads to some mismatch in the FD and the PD mo- 
mentum ranges and this in turn results in a limitation on 
the pp — » ppu! phase space available within the ANKE 
acceptance. The acceptance corrections can therefore not 
be done in a model independent way and some assump- 
tions have to be made on the differential distributions to 
determine the total cross section. To test the validity of 
these assumptions, we have compared a number of one- 
dimensional pp ppco experimental distributions with 
the results of simulations. The following observables have 
been considered: 

— The excitation energy in the two-proton rest frame: 
Epp = Afinv (p,_p) - 2mp. 

— The excitation energy of the w-meson and the proton 
detected in the PD: e^^p = Minv(w,p) — mp — rriu;. 

In this case, the total range of the Dalitz plot pro- 
jection could be covered due to the large FD momen- 
tum acceptance. The range of the other projection, 
where the proton was detected in the PD, was less 
than 40 MeV and 75 MeV at excess energies of 60 MeV 
and 92 MeV, respectively. 

— The angle of the w-meson with respect to the beam di- 
rection in the overall centre-of-mass system: cos^cm- 

— The angle of a proton momentum q with respect to 
the beam direction in the pp rest frame: cos 0pp. 

— The angle of q with respect to the momentum of the 
w-meson in the pp rest frame: cos!?pp. 
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The number of w-mesons in each bin of such a distribu- 
tion was determined from the corresponding missing-mass 
spectrum using approaches that depended sliglitly on the 
shape of background and tlic position of the to peak rela- 
tive to the upper edge of the spectrum. 

— When, as in Fig. 5, the shape of the background could 
be clearly established, a preference was given to the fit 
procedure described above. The missing-mass trans- 
formation was used to check that the description of 
the backgroimd below the lu peaks was reasonable. It 
was found that the normalisation factor of the trans- 
formed missing-mass distribution to the measured one 
is not constant, but may vary by up to 7%. Therefore, 
this factor had to be determined for each bin. 

— In cases where the fit procedure could only be ap- 
plied at a single excess energy, generally at 92MeV, 
the background function for the corresponding trans- 
formed distribution was obtained. A simple scaling of 
this function was then allowed which, together with 
a Gaussian oj peak, was used to fit the distribution 
measured at the other excess energy. The scaling pa- 
rameter thus obtained was used for the normalisation 
of the transformed distribution to the measured one. 
The number of w-mesons obtained as a result of the fit 
was verified using the difference of distributions of the 
type shown in Fig. 4. This procedure was also applied 
at Epp > 60 MeV, where the distributions transformed 
from 2.85GeV/c to 2.95GeV/c no longer contained 
the u) peak. 

— The scenarios outlined above describe more than 80% 
of our data. For the remainder, the transformed dis- 
tribution was normalised by the ratio of the numbers 
determined after excluding the range covered by the 
actual and the transformed a;-peaks from each of the 
spectra. It was verified that the sum of background 
events under the actual w-peak found in this and the 
previous two cases is equal within errors to the total 
number of background events shown in Fig. 4. 

In all the cases, difii'erent assumptions regarding the 
background functions and/or different normalisation fac- 
tors were tested in order to establish possible systematic 
errors, and these are incorporated into the vertical error 
bars on the points presented in Fig. 6. 

The pp ppuj reaction was simulated at both ex- 
cess energies using the PLUTO event gcinerator [14] under 
different assumptions on the distribution of events over 
the total phase space of the reaction. Both protons were 
traced through the ANKE acceptance with the help of the 
GEANT3 package [15]. To check the tracing procedure, 
the momentum and angular acceptances of all the scintil- 
lators, as well as their images in the MWPCs, have been 
evaluated and compared with the corresponding experi- 
mental ones. They were found to coincide to within the 
experimental errors. Finally, the total number of accepted 
events simulated at each excess energy was normalised to 
the corresponding total number of w-mesons obtained in 
the experiment. 

The solid lines in the Fig. 6 demonstrate the distri- 
butions expected in the case of a three-body phase space 
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Fig. 6. The distributions of experimental (points) and simu- 
lated (lines) w-meson yields for different reaction observables 
at 60 MeV (left panels) and 92 MeV (right panels). The differ- 
ent simulations correspond to pure phase space PS (solid line), 
full FSI (dashed line), phase space distorted by an angular de- 
pendence in the variable (dash-dotted line), and partial 
FSI with the distorted angular distribution (dotted line), as 
described in the text. The corresponding total cross sections 
are given in Table 2. 



(PS), whereas those shown by dashed lines take into ac- 
count the final state interaction (FSI) between the two 
protons. Following Goldberger and Watson [16], the FSI 
enhancement factor 7(g) was calculated using the Jost 
function J(g): 



0" 



(4.1) 
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where q = \q\ is the magnitude of the proton momentum 
in the pp rest frame. On the basis of the pp scattering 
Icngtli and effective range we took a = —20.5 MeV/c, (3 = 
167MeV/c. 




Fig. 7. Missing-mass distribution of the pp ppX reac- 
tions obtained at Q = 92 MeV in the 10 < £pp < 30 MeV 
range (closed circles) is shown together with the corresponding 
transformed distribution (crosses). The thick solid line repre- 
sents the fit result for the measured spectrum. The dashed and 
dashed-dotted lines show two extreme cases of the background 
behaviour. The maximum systematic error was established to 
be 25%. 

Due to the selection of one proton in the FD and the 
other in the PD, the experiment had no acceptance for 
Epp < 10 MeV. As a consequence, and as shown by the 
left panels of Fig. 6, the distributions at Q = 60 MeV 
cannot distinguish between the pure PS and FSI mod- 
els. However, at Q = 92 MeV, the experimental yield in 
the 10 < Epp < 30 MeV range is about a factor of two 
lower than that expected on the basis of the FSI model. 
It must be noted that the yield in this bin can be de- 
termined very reliably, as shown by the corresponding 
missing-mass spectrum in Fig. 7. The fit procedure ap- 
plied to the missing-mass distribution at Q = 92 MeV 
would give an even smaller yield than that found taking 
into account the behavior of the transformed spectrum un- 
der the peak. Different background functions were tested, 
but none of them could increase the yield by more than 
25% assuming a smooth shape under the lu peak. 

The low experimental yield at Q = 92 McV for Spp < 
30 MeV is similar to that observed in the analogous pp — > 
ppr] reaction at 74 MeV, where the effects of the FSI in 
the pp system were found to be very small [17]. This 
may be an indication that the pp S-wave is weak and 
that higher pp partial waves are important. In this con- 
nection it should be noted that the distributions in the 
proton angles in the pp rest frame, especially, ^pp, would 
allow some anisotropy. This contrasts with the yield distri- 
butions in the meson centre-of-mass angle, which corre- 



Table 2. Total cross sections for different models. The form 
of the FSI factor is given in eq. (4.1) and that of the assumed 
angular dependence in eq. (4.2). 





Q = 60MeV 


Q = 92 MeV 


Model 


a = 0.2 


a = 0.6 




k = 0.8 


A; = 0.4 


PS 


5.0 ph 


10.8 fih 


PS + FSI(100%) 


7.1 ^b 


14.3 fih 


PS + a P2 (cos "Z'pp) 


4.8 ^b 


10.4 nh 


PS + kxFSl + aP2 (cos %p) 


6.6 iA> 


12.0 /ib 



spond quite well to an isotropic da/dcosOcm at both the 
excess energies. The TOF collaboration reported a similar 
behaviour [8]. 

The dash dotted lines in Fig. 6 were calculated as- 
suming that the three-body phase space is modified by a 
factor 

1 -haPa (cos tf'pp), (4.2) 

where P2(cos!2'pp) is a Legendre polynomial, and a = 0.6 
at Q = 92 MeV. It is interesting that this assumption gives 
a better description of the data, not only for the cos'I'pp 
distribution, but also for that of cos ©cm in the —0.9 re- 
gion. At Q = 60 MeV we would expect any anisotropy to 
be smaller. The case of a = 0.2 was tested but, for all 
the observables other than costf'pp, the results differ little 
from the PS case. 

Although the introduction of an angular dependence 
through eq. (4.2) provides a quite adequate description of 
all the distributions, with no need for the S wave FSI, it 
seems unlikely that this would be suppressed completely. 
To test some of the possibilities, we have adjusted the con- 
tribution of the FSI to provide a slightly better description 
of the data; the results shown by dotted lines in Fig. 6 were 
achieved by decreasing the FSI enhancement by factors of 
0.8 and 0.4 at Q = 60 MeV and 92 MeV, respectively. The 
values of the total cross sections obtained on the basis of 
these various assumptions are presented in Table 2. Other 
variants, including allowing some dependence on 0pp, in- 
variably give results that lie within the ranges shown in 
the table. 
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5 Results and Discussion 

Combining the results of all the variants that can hardly 
be distinguished within the errors, we find total cross sec- 
tions lying in the ranges 4.8 — 7.1/ib at Q = 60MeV and 
10.4 - 14.3^b at 92McV, respectively The highest values 
are achieved with the rather unrealistic pure FSI model 
and so we take rather the last two lines of Table 2 as an 
indication of the model-dependence of the results. These 
give 

a{Q = 60 MeV) = (5.7 ± O.e^tat ± O.Ssyst ± O.g^od) Mb, 
a{Q - 92 MoV) = (11.2 ± l.lsiat ± 1.7syst ± 0.8„,od) ^ih , 
where the first error represents the statistical uncertainty. 
The second figure stands for the total experimental sys- 
tematic imccrtainty, which includes contributions from the 
luminosity determination, the possible change of accep- 
tance due to uncertainties in the detector positioning, 
but mainly the determination of the total number of ppuj 
events. These were all treated as independent errors. The 
final number is our estimate of the influence of the model 
on the acceptance calculation. 




Q [MeV] 



Fig. 8. The energy dependence of the pp ppui total cross 
section. Our two points (closed circles) arc compared with ear- 
lier data from Refs. [6] (triangles), [7] (square), and [8] (open 
circles). In all cases, systematic and statistical eros have been 
added in quadrature. The chain and dashed curves show re- 
spectively the normalised energy dependence of eq. (5.4) with 
and without smearing over the finite ui width. The solid line is 
a purely phenomenological interpolation of the whole data set 
given by eq. (5.3). 



where Q is measured in MeV and the cross section in jjh. 
This is also shown in the figure. 

In the SPESIII work [6], the energy dependence of the 
total cross section was compared with a simplistic model, 
consisting only of three-body phase space but modified by 
the 5-wave pp FSI [18]. This leads to the analytic result: 

.ta-^) = a (^^-^J, (5.4) 

where e « 0.45 MeV. However, at low excess energies it 
is important to smear this function over the width of the 
w-meson. Both these forms arc shown in Fig. 8, with the 
value of — 37 nb being chosen to fit the SPESIII points 
rather than those at higher energy. Although the ^-wave 
assumption is clearly untenable at larger Q, the curves do 
suggest that P or higher partial waves are important there 
and a good representation of the data can be obtained by 
adding a small amount of pure phase-space distribution 
to the FSI of eq. (5.4). While presenting a baseline against 
which one can judge the more refined theoretical models, 
this approach also shows that the finite-width effects are 
negligible in our energy domain. 

The predictions of three microscopic models [3-5] are 
shown in Fig. 9. All of them are versions of a meson- 
exchange model, but with different emphasis on the role of 
the N* isobars. It should be noted that these calculations 
involve a large number of parameters, some of which have 
been adjusted in order to optimise agreement with the 
published data, including the differential distribution at 
173 MeV [8]. 
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Our pfp — > ppio total cross sections are shown in Fig. 8 
with vertical error bars estimated from the quadratic sum 
of both the systematic and statistical errors presented 
above. Also shown are the results of other experiments at 
low Q [6-8]. The 92 MeV value is slightly higher than that 
of TOF, though the error bars do overlap. With the possi- 
ble exception of the highest SPESIII point, the results are 
consistent with a smooth energy dependence and a global 
fit to the whole data set yields 

logio(T= -0.855 + 0.495 logiog + 0.230 (logioQ)', (5-3) 



Fig. 9. Experimental data on the pp — > ppw total cross section 
shown in Fig. 8 compared with the predictions from micro- 
scopic theoretical models shown by the chain ([3]), solid ([4]), 
and dashed ([5]) curves. 

In the early work of Tsushima and Nakayama [3] , nu- 
cleonic and mesonic current contributions were considered 
but the energy dependence of the total cross section is de- 
scribed better by the inclusion of effects of nucleon reso- 
nances, and it is that curve which is shown in Fig. 9. Nev- 
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ertheless, this calculation still underestimates the SPE- 
SIII data by about a factor of two even though the finite 
LO width was taken into account. Despite neglecting this, 
and the initial state interaction, Kaptari and Kampfcr [4] 
got good agreement with the then available data without 
relying on any N* contribution. It was shown there that 
the effects of the 5'-wave FSI are particularly important, 
for without it a phase-space energy dependence of the 
type is predicted. The Tiibingen group [5] had exten- 
sive recourse to several nucleon resonances, especially the 
A''*(1535), but the novel feature of their approach is the 
coupling to off shell lu mesons, which is particiilarly im- 
portant close to the reaction threshold. If this is correct 
then one would expect interesting effects when the oj is 
measured exclusively through its e+e^ decay mode [5]. 

Having established the energy dependence of the total 
cross section for pp — > ppu), we are now in a position to 
investigate further the OZI rule in near-threshold produc- 
tion processes. If the (w, (j)) mixing in the vector meson 
nonet were ideal, such that the (j) were a pure ss state, 
then its three-pion decay would be forbidden by the OZI 
rule [2]. That it takes place at all is interpreted as an 
indication that the mixing is not quite ideal. Using the 
Gell-Mann-Okubo mass formulae, Lipkin [19] estimated 
that the ratio of coupling strengths of the ^ and ui to 
non-strange hadrons should be 



i?ozi = 4.2 X 10" 



(5.5) 



It is therefore useful to evaluate the ratio of the (j) and ui 
production cross sections in pp collisions at the same value 
of the excess energy: 



R<t>/u> = 



_ a{pp pp4>) 
(j{pp PP^ ) 



(5.6) 



Taking the recent measurement [9] of the total cross 
section for (j) production in pp collisions at Q = 18.5, 
34.5, and 75.9 MeV in conjunction with the global fit of 
eq. (5.3) to the data set of Fig. 8, we find at these ener- 
gies R^/^ = (3.1 ± 0.6) X 10-2, (3.0 ± 0.7) x 10"^, and 
(2.4 ± 0.7) X 10"^, respectively, where all error bars have 
been compounded quadratically. Applying the same pro- 
cedure to the (j) point reported by the DISTO collabora- 
tion leads to R^/^ = (2.2 ± 1.0) x 10"^ at 83 MeV. The 
mean value of -R^/^^ = (2.8 ± 0.4) x 10"^ sa 7 x i?ozi, is 
lower than that obtained in Ref. [9] because of the higher 
values for lj production that are now apparent in Fig. 8. 
Although not statistically significant, there are indications 
from these numbers that R^ might decrease with rising 
Q but; to test this, data would be needed on (j) produc- 
tion at higher values of Q. However, it should be noted 
that such a slow decline had been predicted in model cal- 
culations. This arises there through the decreasing effect 
of the destructive interference between the nucleonic and 
mesonic contributions to u production as the energy is 
raised [20]. 



6 Conclusions 



We have presented new measurements of the pp ppui re- 
action at excess energies of 60 and 92 MeV. The SPESIII 
technique [6] of kinematically shifting the data at one en- 
ergy to estimate the background at a neighbouring energy 
was successfully employed for differential as well as total 
spectra. In this way the small u) signal in the missing-mass 
distribution could be identified despite the large amount of 
multipion production. With the setup actually employed 
for studying this reaction at ANKE, the acceptance, espe- 
cially at small cpp, was limited. A variety of assumptions 
were tested against one-dimensional spectra in order to 
assess the model dependence of the total lo production 
cross section that we obtained in this way. This led to 
uncertainties that were comparable to the statistical and 
other systematic errors. The value for the total cross sec- 
tion at 92 MeV is a little higher than that found by the 
TOF group [8], but our two points join smoothly with the 
results of other experiments. There is evidence from the 
anisotropic pp angular distribution at 92 MeV that P or 
higher pp partial waves are significant at this energy. 

The new data allow us to extract the OZI ratio with 
greater confidence and smaller error bars than before. This 

results in a slightly smaller value of than previously 
quoted [9] , with possibly a hint of some energy dependence 
(but see also Ref. [21]). 

All three theoretical calculations considered [3 5] ap- 
peared after the broad outlines of the energy dependence 
of the total cross section had been determined experimen- 
tally and in some cases the model parameters were ad- 
justed to reproduce this. Other experimental observables 
such as differential cross sections or decay angular distri- 
butions are needed to constrain the models more tightly 
and it is encouraging to note that analysing power results 
will soon be made available by the COSY-TOF collabo- 
ration [22] . Since good data exist on pn — > d(j) at low en- 
ergies [23], further data on uj production with a neutron 
target would be particularly helpful [24] . 

The WASA at COSY facility will soon become opera- 
tional [25]. The production of the u) in pp collisions could 

then be investigated, with larger geometrical acceptance 
and lower background, via the detection of the lo — ^ n^j 
decay {BR w 9%) and the reconstruction of its invariant 
mass. Such an approach would also have the advantage 
of leading to a determination of the tensor polarisation of 
the LO, a quantity which is very sensitive to the presence 
of higher partial waves. A similar sensitivity also exists 
in the spin-spin correlation in P ppLO, which could be 
measured at ANKE [26]. 
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